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The dynamics of self-oscillatory extended systems, resonantly forced at a frequency close to 
that of the natural oscillations (1:1 resonance), is shown to be universally described by a complex 
Ginzburg-Landau equation containing an inhomogeneous term. The case of amplitude modulated 
forcing is considered, which generalizes previous studies. Application to the two-frequency forcing 
in a 1:1 resonance is considered as an example. 
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O ■ I- INTRODUCTION 

(N , 

The periodic forcing of spatially extended self-oscillatory systems is a classical method to excite the formation of 

, spatial patterns in such systems. This kind of forcing admits a universal description when the system is operated near 

I ■ the oscillation threshold, and forcing acts on a n : m resonance, defined by the relation oje — (n/m) (luq + i^) between 

I the external forcing frequency tUe and the natural frequency of oscillations ujq, where n/m is an irreducible integer 

I— I . fraction and i/ is a small mistuning. In such case the slowly varying complex amplitude of the oscillations u verifies 

g5 . the following g„e.ali„d co„pl„ Gi..b„g-L„da„ e,„at.o. (CGLE, § 

' dtu = aiu + a-zV'^u + a3\u\ u + a4u"' ^, (1) 

where u stands for the complex conjugate of u, 0^=1^2,3 are complex coefficients and 04 is proportional to the mth 
power of the forcing amplitude. Eq. M) is valid in principle for perfectly periodic forcings and has been introduced 
making use of symmetry arguments . That elegant reasoning is based on the existence of a discrete time invariance 
' when a perfectly periodic forcing is applied to a system. However if, e^., the amplitude of forcing is modulated in 
' time in an arbitrary way all temporal symmetries are broken and Eq. (|l|) is not rigorously justified. In this work we 
prove that Eq. does provide a universal description of the close to threshold dynamics of self-oscillatory extended 
' systems forced in a 1:1 resonance {n = 1 and m = 1 in Eq. (|l])). The analysis is based on the technique of multiple 
scales and generalizes the concept of resonant forcing as it considers (almost) periodic forcings that are nonuniform 
across the system and/or modulated in time. Our derivation hence allows a rigorous simplified study of both noisy 
forcings as well as multi- frequency forcings within a 1:1 resonance. Generalizations to other resonances (n ^ 1 or 
m ^ 1), although cumbersome, can be made straightforwardly following the lines of the present derivation. 
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II. MODEL 

We consider a generic two-dimensional system described by N real dynamical variables {Ui (r,i)}^i whose time 
evolution is governed by the following set of real equations written in vector form, 

dtV{r,t) = i{fi,a;U,V^V), (2) 

where f is a sufficiently differentiable function of its arguments. We assume that the dependence of f on spatial 
derivatives of U is through its Laplacian = + dy, r = [x,y). This is the simplest dependence on derivatives 
in rotationally invariant systems at the time it corresponds to physical systems of most relevance like, e.g., reaction- 
diffusion and nonlinear optical systems, /i is the bifurcation parameter and a (r, t) is the forcing parameter, which is 
allowed to vary on time and space. Physically a may represent either an independent parameter, or the modulated 
part of any other parameter. 

We assume that in the absence of forcing {a = 0) Eq. (||) supports a steady, spatially homogeneous state U = Us (/i) 
((9fUs — dx^s — dylJs = 0), which looses stability a,t — fio giving rise to a self-oscillatory, spatially homogeneous 
state. In other words, we assume that the reference state Us suffers a homogeneous Hopf bifurcation dX ^ — hq. 
We wish to study the small amplitude solutions that form in the system close to the bifurcation when the arbitrary 
parameter a is modulated in time with a frequency close to that of the free oscillations. 

For the sake of convenience we introduce a new vector 
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u(r,t)=U(r,i)-U„ (3) 

which measures the deviation of the system from the reference state, in terms of which we rewrite Eq. (|^) as a Taylor 
series, 

dtu ^F{n,a)+J{n,a)-u + V (^, a) ■ V^u 

+K {fi, a- u, u) + L (/i, a; u, u, u) + h.o.t., (4) 

where h.o.t. denotes terms of higher order than 3 in u or than 1 in V^u. These h.o.t. have no influence near the 
bifurcation (/x /ig) since, as we show below, they are O — /io|^^ or smaller and only terms up to O ~ Mol'^^^^ 
contribute to the leading order dynamics of the system whenever the bifurcation is supercritical, which is the case we 
assume. 

The different elements of the expansion (^) are defined as 

F(^,a) is, 

ifi, a) - [df,/dU,]^ , A, ifi, a) = [df,/dV'U,] ^ , 

K ifi, a- a, b) = if E5=i [d^i/dU,dUj] ^ a,b„ 

L(^,a;a,b,c) = ^Yl,tj,k=i [d^^ / dUidU^dUk] ^a.bjCk, 

where a, b,c are arbitrary vectors and the subscript s denotes U = Us (/^). Vector F (/i, a) is subjected to the 
condition 

F(m,0)-0, (5) 

since in the absence of forcing [a — 0) the reference state u = is a steady state of Eq. by hypothesis. J and V 
are matrices, and vector K (L) is a symmetric and bilinear (symmetric and trilinear) function of its two (three) last 
arguments. 

III. THE HOPF BIFURCATION 

In the absence of forcing the stability of the reference state against small perturbations 5u is governed by the 
following equation 

^tJu = J (^, 0) • (5u + P (^, 0) • V^^u, (6) 

obtained upon linearizing Eq. (|^) for a = with respect to 5u. The general solution to Eq. (^ is a superposition of 
plane waves of the form 



(r, t) = Y.j Wj exp {K.jt) exp (ik^-r) , (7) 



with 



A,w, =7U(^,A:|) -w,, (8) 
M {li, k") = J ifi, 0) - k^V (^, 0) , (9) 

{k"^ = k • k) hence eigenvalues and eigenvectors of matrix M depend on k only through fc^ as A4 does. As we are 
assuming that the reference state looses stability at fi = fiQ via a homogeneous Hopf bifurcation in the absence of 
forcing, matrix M. {fJ-jk"^) must have a pair of complex-conjugate eigenvalues {Ai, A2} — {A [^J-,k^) , A (the 
overbar denotes complex conjugation) governing the instability, i.e.: 

(i) Close to the bifurcation ReAi>3 < whilst Re A can become positive for some fc's, 

(ii) At the bifurcation Re A is maximum and null at fc = (the perturbation with largest growth rate is spatially 
homogeneous): 

ReAo = 0, (afcReA)o = 0, (S^ReA)^ < 0. (10) 

where, here and in the following. 
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subscript affecting functions denotes {fi = iiq, a = 0, k = 0} . 

(iii) The instability is oscillatory, i.e., 

lmXo=ujo^O. (11) 

Finally, the preceding properties imply that: 

(iv) All eigenvalues of Mo ~ Jo, see Eq. (^), have negative real part but {Aq, Ao} = {it^o, —iuJo}- 

For the sake of later use we introduce the right and left eigenvectors of JTq associated with eigenvalues {iluq, —iluq}, 

ht. Jo = «c^oht, ht. J'o = -*woht, ^ ' 

where the short-hand notation h = wi (/i = /io, A;^ = O) , h = W2 (/i = /io, fc^ = O) has been introduced. These vectors 
verify the following orthonormality relations: 

• h = 0,h^ ■ h = 1. 



IV. SCALES 



We are interested in determining the small amplitude solutions that from in the system close to the Hopf bifurcation, 
which we define by 



(13) 



where e is a smallness parameter (0 < £ <C 1). The study is based on the widely used technique of multiple scales. 
These spatial and temporal scales appear naturally close to the bifurcation and are those on which the asymptotic 
dynamics of the system evolves. As is well known, in a homogeneous Hopf bifurcation these slow scales are given by 



T^eH,X = sx,Y = ey. 



(14) 



These scales follow from the behaviour of A close to the bifurcation, Eq. ( |13D , for values of k close to the most unstable 
mode k — 0: 



1 



A (/io + £^2, A:') = Ao + (a^A)o £^2 + ^ {dlX}^ ^ + max {O {e^) , O {k") ] 



(15) 



where the term {dkX)^ k has not been included since ((?fcA)Q = as A is an even function of k. From this equation we 
obtain, making use of Eq. (Eol), 



Y^eX{^lo + e^^l2,k^) = (a^ReA)o£2/i2 - - Re A|q fc2 + max {£)(£'*) ,0(fc2)} , 

what indicates that the only modes which can experience linear growth verify k — O {e), hence the asymptotic dynam- 
ics of the system exhibits spatial variations on a scale x, y ~ k^^ ~ £^^and the slow spatial scales {X ~ ex, Y = ey) 
follow. Thus, putting k — eki. 



ReA (/io +£^M2,e^fci) = £^ 
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(a^ReA)„A.2-:.|5,^ReA|„fcf 



which shows that the growth of perturbations occurs on a scale t (Re A) ^ ^ £ ^ and the slow timescale T = e^t 
follows. On the other hand, making use of Eq. (pT|) and putting k — ek\ again, Eq. (fst) yields 



Im A (/io -!-£'' /i2,e^fci'' 



(a^ImA)„/i2 + -(a2lmA)„fc2 



0\ 
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whose first term, ujq = O (^e^) , indicates that the original timescale t must be retained, whilst the rest of terms do not 
introduce other relevant timescales. 

As for the external forcing we assume that its form is consistent with the previous scales. In particular we assume 
that it is weak and of the form 

a(r,t) = £^^3 (X, F, T, t) , (16) 

^ [aa.p {X, Y, T) exp {iptoot) + ag.p (X, F, T) exp {-ipuj^t)] , (17) 
p— 1 

Thus we are considering a resonant forcing with slowly varying amplitude. Notice that, according to Eqs. ( p^pT| ), 

a l^r, i + = a (r, t) + O (e^) and hence the considered forcing is almost periodic of fundamental frequency Wg- 

Under these conditions a multiple scale analysis is possible and we look for asymptotic solutions to Eq. (|^) in the 
form 

u(r,i)=V°° e"u™(X,r,r). (18) 
We finally introduce Eqs. and (|l6|-[l8|) into Eq. ^ making use of the following chain rules for differentiation 

^ — ^m— 1 ^ ' 
^ — ^m— 1 ^ 

and solve at increasing orders in e. 



V. THE COMPLEX GINZBURG-LANDAU EQUATION 

The general form of Eq. (^ at any order e™ is found to be 

3(u™) = g™(X,y,T,i), (19) 

where 

a (u) = dtu ~Jo-u, (20) 

and gm does not depend on u™ (but on u„<m). Clearly, as 5 (exp {ituot) h) — ^ (exp (—iujot) h) — 0, see Eq. (|T^), 
the solvability of Eq. ( p^ ) requires 

dt'h^ ■ g™ (X, Y, T, t') exp (-^wo^') = 0, (21) 
(or its equivalent complex-conjugate) which ensures that g^ does not contain secular terms (proportional to 



exp {iuJot) h or to exp {—iujat) h). Once condition (|2l| ) is verified, the asymptotic solution to Eq. (19) reads 

Urn [X, r, T, t) = Urn {X , Y, T) exp {iLOot) h + u„ (X, F, T) exp {-iujot) h 

+u4(X,y,T,t), (22) 

where Um {X, Y, T) is a function of the slow scales, and the last term is the particular solution. Note that the solution 
( p^ ) should involve, in principle, terms proportional to all the eigenvectors of Z (•) [which are those of Jo, see Eq. (||)]. 
However all of them are damped according to exp [— |Re (/^o, 0)| i], since Re Ai>3 = /iq, A: = 0) < by hypothesis, 
except those associated with (h, h) . 

A. Order e 

This is the first nontrivial order and reads 

gi = 0. (23) 



The solvability condition (21) at this order is fulfilled and, according to Eq. (|22|), 



ui = Ml (X, y, T) exp {iujot) h + ui (AT, F, T) exp {-iuoot) h. (24) 
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B. Order 



At this order 



g2 (5pF)q^2 +K(/io,0;ui,ui). 



(25) 



The first term of the r.h.s. is null by virtue of Eq. (^). Making use of Eq. (|2^) and taking into account that K is 
symmetric and bilinear in its two last arguments, Eq. (25) can be written as 



g2 = 2K (^0,0; h,h) |Mi|' + K(/io,0;h,h)M2exp(i2woO 
+K (^0, 0; h, h) u\ exp {-i2ujQt) . 



The solvability condition (21) is automatically fulfilled again and, according to Eq. (|22|), 

U2 U2 {X, Y, T) exp {iuot) h. + u2 {X, Y, T) exp {-iuot) h 



where, 



-vo |ui I + V2U^ exp {2iujQt) + V2U^ exp {—2iujQt) . 

Vo = -2Jo-i-K(Ato,0;h,h), 

V2 = - (Jo - i2uJoI)~^ ■ K (mo, 0; h, h) , 



(26) 

(27) 

(28) 
(29) 



are constant vectors, and T is the N x N identity matrix. Note that both Jo and {Jo — i^ojol) are invertible since 
neither nor 2iujo are eigenvalues of Jo by hypothesis: otherwise other eigenvalues different from {iwq, —iujo] would 
have null real part at the bifurcation. 



C. Order 



Finally, at this order we find 



g3 = -drui + ^2 (5^ J)o • ui + • (5i + af.) ui 

+ (,daF)o as + 2K (/xq, 0; Ui, U2) + L (^0, 0; Ui, Ui, Ui) . (30) 



Application of the solvability condition ( |2l| ) yields, after substituting Eqs. (17), (|2J) and ( pTj ) into Eq. (|30|), and 
making use of the symmetry and linearity properties of vectors K and L, 



where 



dxui = C1H2U1 + C2 {d\ + dy) ui + C3 ui + C4a3,i, (31) 

ci = ht . (d^J)^ ■ h, (32) 

C2 = . Do • h, (33) 

C3 = 2ht . K (^0,0; h, vq) + 2ht • K (/xq, 0; h, V2) + 3ht • L (/zq, 0; h, h, h) , (34) 

C4 - ht . (a„F)o , (35) 

are constant coefficients. 

Finally note that, making use of Eqs. (p^), ( p^ and (|2^), the asymptotic state of the system can be written as 

u (r, t) =u (r, t) exp (icjoO h + u (r, t) exp {-iuJot) h + O (e^) , (36) 
u{r,t) =eui{X,Y,T,t), 

hence u denotes the leading order amplitude of oscillations. Its evolution equation is obtained from Eq. (51^) by 



returning to original scales and parameters via Eqs. ( |13| ) and (14), and reads 

dtu (r, t) = (ji- ^0) cm + C2V^u + C3 \uf' u + f {r,t), (37) 

where 
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fiv,t) = c^e''a3,i{X,Y,T) 



(38) 



is a slowly varying function proportional to the complex amplitude of the fundamental component of forcing e^aa^i, 
see Eqs. (16|, |l7[). 

Eq. (p7), orTsy), is a complex Ginzburg- Landau equation containing an inhomogeneous forcing term, / (r, t), which 
generalizes Eq. (0) for n = m ^ 1. Eq. (|3^ ) is valid whenever Recs < (supercritical bifurcation) since otherwise 
it could lead to unbounded solutions. If Re C3 > the bifurcation is subcritical and the analysis must incorporate 
higher orders in the e-expansion. 



VI. APPLICATION TO TWO-FREQUENCY FORCING 

Just for the sake of illustration let us finally consider the two-frequency forcing case. In particular we assume that 
the forcing parameter a has the form 

a{T,t) = A{u;it) + A{uj2t), (39) 

where ^ is a 27r periodic function [A {0) = A{0 + Stt)]. This means that we are dealing with the superposition of 
two forcings of equal amplitude but of different frequencies (wi and ^2)- An example of this type of forcing could be 
the illumination of a photosensitive version of the Belousov-Zhabotinsky reaction with two equal sources, whose light 
intensities are periodically modulated in time at two different frequencies. Another example could be the injection 
into a laser cavity of two coherent fields of equal amplitudes and different frequencies. 
Due to the commented periodicity one can write 

A [ujit] ~ ^ [-^p 6xp (ipijJit) + Ap exp {~ipijjit)\ . (40) 

(Note that the term p = is absent, as in Eqs. ( [l6| , |l7|) , since it would correspond to a constant bias. This bias term, 
if any, is implicitly considered in our analysis at the same level as any other constant parameter of the system and it 
would appear, in principle, in the determination of the Hopf bifurcation). Upon rewriting Eq. (^) as 

A {ujit) = ^ _ [Ap exp {ip6it) exp (ipuJot) + Ap exp {—ip5it) exp {-ipLO^t)] , (41) 
where 5i — uji — ljq^ we can express Eq. ( |3^ in the form ( |l^Jl7| ): 

E°° 
[ap (t) exp {ipujot) + ap {t) exp {-ipujot)] , (42) 
p—i 

where 

ap (t) — Ap [exp (ipSit) + exp {ip62t)] — 2Ap exp (ipi/t) cos {ptot) , (43) 

where 

di + S2 LDl + UJ2 Si — S2 UJl — W2 

-=^^ = ^—--0, 0.= ^— = ^^. (44) 

(We note incidentally that, in the case of single-frequency forcing, uji ~ lu2 = uje and hence v = oj^ — liJo {oJe =^^0 + ^^)1 
and Lo — 0.) Hence all our previous analysis is valid for this type of driving, whenever forcing is weak {Ap ^ e^) and 
the time scale along which the a^'s vary is slow. Specifically this requires i^, w ~ e^. If these relations are satisfied, 
the system will be described by Eq. ( ^ with / given by Eq. ( ^ and s^a^s given by ap=i {t): 

dfU (r, <) — (/i — iiq) ciu + C2^^u -\- C3 \u^ u + 2C4A1 exp {ivt) cos {uit) . (45) 
In order to simplify this equation let us define 

C/(r,i) = M(r,i)exp(-ii/i) , (46) 

which, substituted into Eq. ( p5[ ) yields 

dtU (r, t) = - ^io) ci -ii^]U + czV^fJ + C3 |?7|^ U + 2ciAi cos (ut) . (47) 
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The actual state of the system, Eq. (|36|), is given in terms of U as 

u (r, t) = 2 Re {[/ (r, t) exp [i (ujq + i') A h} . (48) 

Finally, note that, especially in nonlinear optics, one is used to express the state of the system in terms of the 
negative-frequency part of the oscillations, i.e., as 

u (r, i) = 2 Re { [7 (r, t) exp [-i {loq + v) t] h} , (49) 

and the complex amplitude U (r, t) verifies the complex conjugate of Eq. (^): 

dttJ{v,t) = - iiq)ci + w]U + C2V^U c:i\U\^ U ^ 2ciAi cos {ujt) . (50) 

This absolutely trivial transformation is done here just to point out that the sign that affects the mistuning v in the 
forced CGL equation depends on whether one uses the amplitude of the positive- or the negative-frequency part of 
the oscillations in order to describe the dynamics of the system. 
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